FTUV-11-0620 KA-TP-13-2011 LPN11-32 SFB/CPP-11-31 
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In this letter we report on a calculation of W 77+jet production at next-to-leading order QCD. 
We include the leptonic decays of the W and take into account all off-shell and finite width ef- 
fects. This is the first computation which falls into the category of triboson+jet production at 
next-to-leading order QCD. In total we find sizable corrections with nontrivial phase space depen- 
dencies. Therefore, our results are important for phenomenological analyses such as the extraction 
of anomalous electroweak quartic couplings from inclusive hadron collider data. 



I. INTRODUCTION 



The production of multiple electroweak bosons is an 
important channel to test experimental data, collected 
at both the Tevatron and the Large Hadron Collider 
(LHC), against the theoretically well-established Stan- 
dard Model (SM) hypothesis. As the mechanism of elec- 
troweak symmetry breaking is yet to be determined, pre- 
cise predictions of electroweak gauge boson production 
rates are necessary to experimentally infer deviations 
from the electroweak coupling pattern predicted by the 
SM. It is well-known that computations of production 
cross sections and differential distributions suffer from 
severe theoretical shortcomings if they are limited to the 
semi-classical (i.e. leading order, LO) approximation in 
perturbation theory. The arising uncertainties intrinsic 
to fixed order calculations are conventionally assessed by 
investigating variations of renormalization and factoriza- 
tion scales, which are remnants of the perturbative series' 
truncation. 

The next-to-leading order (NLO) real emission contri- 
bution to the hadronic cross section, however, can probe 
new partonic initial states not present at LO, so that the 
LO scale uncertainty can sometimes be totally mislead- 
ing. This is especially true for processes which are char- 
acterized by a QCD singlet final state at LO, e.g., elec- 
troweak triboson production [TJ [2J. For these channels 
the total NLO correction factors are known to be partic- 
ularly sizable, K = er NLO /<7 LO ~ 1.8. The main reason 
for this large correction is that the NLO corrections in- 
clude new, large LO real emission subprocesses initiated 
by gluons. Probing the protons' gluon parton distribu- 
tion at small momentum fractions with the real emis- 
sion contribution, the large correction does not signal a 
breakdown of perturbation theory, but strongly asks for 
perturbative improvements of the one-jet-inclusive cross 
sections as a major contribution to the full next-to-next- 
to-leading order cross sections. Similar observations and 



conclusions hold for diboson [3j 0] and diboson+jet pro- 
duction [SHE] (see also Ref. [§] for a related discussion of 
Z+jets). 

In this paper we report on the first calculation per- 
formed in the context of triboson + jet production: 

PP ^ t'vai + jet + X and P P ^ £+vai + jet + X. 
We include all off-shell and finite width effects in our 
calculation, i.e. we compute the full matrix element at 
0(a ot 2 s ). For convenience we will refer to these processes 
as W 77 + jet production. 

We organize this paper as follows: Section [TT] reviews 
the technical details of the calculation and comments on 
the numerical Monte Carlo implementation. In Sec. |III| 
we discuss the numerical results; we examine the cross 
sections' scale variations and the impact of the QCD cor- 
rections on differential distributions. Section [Tv| gives a 
summary of this work. 



II. ELEMENTS OF THE CALCULATION 

The leading order 0{a A a s ) contribution to, e.g., V P 
£~£^77 + jet is given by the subprocesses 



qQ -> £~Pe71 + g, 
gQ > £~^77 + g, 
qg -t l~ £^77 + Q , 



(la) 
(lb) 
(1c) 



where q = (d, s) and Q = (u, c) denote the light down- 
and up-type quark flavors, respectively. Due to unitarity 
of the CKM matrix, any dependence on the CKM matrix 
elements drops out for the flavor-summed gluon-induced 
subprocesses of Eq. (lb), (lc). These subprocesses also 
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dominate the hadronic cross sections at the LHC because 
the protons are typically probed at small momentum 
fractions for inclusive production. Consequently, a non- 
diagonal CKM matrix decreases our LO result only at the 
per mill-level. This is well below the residual (NLO) scale 
uncertainty and we therefore use a diagonal CKM matrix 
in our calculation. Furthermore, we do not include bot- 
tom contributions to the hadronic VF77+jet production 
cross section. They either involve top or bottom quarks 
in the final state and are distinguishable experimentally 
by b tagging. 
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FIG. 1: Selected topologies contributing to VF~77+jet pro- 
duction at NLO; i = 1, 2 denotes the generation index. Note 
that we do not show the contributions W„ — > v e + (n — 2)7, 
2 < n < 4. These lead to identical configurations from the 
QCD point of view. The red loops indicate the virtual con- 
tributions, giving rise to topologies of up to boxes (a,b), up 
to pentagons (c) and up to hexagons (d). 



For the numerical implementation of the LO cross sec- 
tion we use routines that are provided by the Vbfnlo 
package [TU] as the real emission contribution to W77 
production at NLO QCD [2]. The hadronic part of the 
amplitude is based on the spinor helicity formalism of 
Ref. [TT], and the electroweak part of the amplitude is 
provided by a cache system which employs MadGraph- 
generated Helas routines [12l H2] using the technique 
of "leptonic tensors" [Hj . All amplitudes of Eq. ([lj are 
related by crossing symmetry, and we show a sample of 
graphs contributing at NLO in Fig. [l] 

For the virtual contributions we use the routines com- 
puted in Ref. [J5] which employ FeynCalc [Hj and 
FeynArts [T7] in an in-house framework. We ap- 
ply the effective current approach described in, e.g., 
Refs. [21 [6l [El [18] , and combine all QCD corrections to 
sets of topologies with two, three, and four attached (ef- 
fective) gauge boson polarization vectors to standalone 
numerical routines.* The order of the gauge bosons is 
thereby fixed and the full amplitude is obtained by sum- 
ming over all allowed permutations. The effective W* 
polarization vectors, see Fig. [I] encode the finite width 



effects in the fixed- width scheme of Ref. [2T] and the off- 
shell contributions of the full decays W* — > tv + (n — 2)7, 
2 < n < 4, in a straightforward way. To be more 
concrete, Figs, [I] (a) and (b) contribute to a routine 
which includes all corrections up to boxes (i.e. all correc- 
tions to qQ — > Wig with terms proportional to Ca and 
Cp — \Ca), Fig. [ll (c) is part of a routine which includes 
corrections up to pentagons (Qg — > W^jq), and Fig. [l] 
(d) contributes to a routine which also includes hexagon 
diagrams (qQ — > W^jjg). The fermion loop contribu- 
tions are sketched in Fig.[2jand were computed within the 
in-house framework [15j and checked against an indepen- 
dent implementation based on FeynArts, FormCalc 
and LoopTools [JJl [22J [23J . 

For the reduction of the tensor coefficients up to boxes 
we apply the Passarino-Veltman approach of Ref. |24) . 
and for a numerically stable implementation of 5- and 6- 
point coefficients we use the Denner-Dittmaier scheme 
laid out in Ref. |25j with the setup and notation of 
Ref. [T5]. This approach has turned out adequate in 
a series of Feynman graph-based hexagon calculations, 
e.g., in Ref. |26j . For completeness, we note that uni- 
tarity cut-based methods have been demonstrated to be 
highly competitive in Ref. |27) . 

The real emission contribution, P P — > t~i>tf) + 2 jets, 
is based on the existing implementations of Refs. [6l[7]- 
We use the C at ani- Seymour dipole formalism [25] to nu- 
merically regularize the soft and collinear QCD diver- 
gences and include finite contributions after cancelling 
the infrared poles of the virtual matrix element |15j to 
order 0(a 4 o^). The real emission matrix element again 
implements the spinor-helicity formalism of Ref. [TT] and, 
analogous to the LO amplitude, the implementation in- 
cludes a cache system for the electroweak W* currents to 
minimize computation time. The implementation of the 
dipoles is optimized to avoid redundancy, and the finite 





"This approach allows to straightforwardly generalize our com- 
putation to models beyond the Standard Model as described in 
Refs. I7i ri9ll20| . 



FIG. 2: Fermion loop contributions to W77+jet production 
at NLO; i = 1, 2, 3 denotes the generation index, q stands for 
all (anti)quark flavors of the incoming (anti)proton. Note that 
analogous gg*y triangle contributions are forbidden by Furry's 
theorem. Not shown are topologies where the polarization 
vector W2 is attached to the initial state (anti)quark, and 
where the internal and external gluon are attached at opposite 
corners of the box. 



3 




FIG. 3: Scale variation of the £ ^77+jet production cross sections at the LHC (£ = e, fi). The cuts are described in the 
text and we take the invariant IV77 mass mw~i~f as central dynamical reference scale. The left panel shows the variation of 
the LO and NLO W /+ 77+jet production cross sections when we change only the factorization scale, only the renormalization 
scale, or both jointly. For W~ 77+jet production the right panel shows the individual contributions to the NLO cross section, 
as discussed in the text. Here we also present results where we have applied a veto on events with two identified jets having 
both a transverse momentum larger than 50 GeV. 



a LO [fb ] 0-NLO jfoj = a NLO /a LO 

W 7 ^ 77+jet 1.191 1.754 1.47 Tevatron 



W+77+jet 4.640 6.634 1.43 
W-77+jet 3.803 5.644 1.48 



LHC 



TABLE I: Total LO and NLO cross sections and K factors 

(-) (-) 
for P V e ^77+jet+X and PP-> e + v e yy+jet+X at the 

Tevatron and at the LHC. The renormalization and factor- 
ization scales are chosen as = /x/? = mw-yy Relative 
statistical and numerical stability errors are below the per 
mill level. 



collinear remainder, which is left after renormalizing the 
parton densities, recycles the born-level matrix elements 
of the dipoles' evaluation and is integrated over the real 
emission phase space applying the phase space mappings 
of Ref. US. 



III. NUMERICAL RESULTS 

We use CT10 parton distributions [53] with a s (mz) = 
0.118 at NLO, and the CTEQ6L1 set [30] with a s {m z ) = 
0.130 at LO. We choose m z = 91.1876 GeV, m w = 
80.398 GeV and G F = 1.16637 x 10" 5 GeV -2 as elec- 
troweak input parameters and derive the electromagnetic 
coupling a and the weak mixing angle from Standard 
Model-tree level relations. The center-of-mass energy is 
fixed to 14 TeV for LHC and 1.96 TeV for Tevatron col- 
lisions, respectively. We consider decays to the two 
light lepton flavors, i.e. for the distributions shown in 
Figs. 3|6 the decays W — > ev e ,[iv^ have been summed, 
and we treat these leptons as massless. 



To study the impact of the QCD corrections on the 
process in detail, we choose very inclusive cuts and a 
strictly isolated photon. A naive isolation criterion for 
the partons and the photon spoils infrared safety by lim- 
iting the soft gluon emissions' phase space. Yet, isolation 
is necessary to avoid non-perturbative jet-fragmentation 
contributions, which would amount to the introduction 
of an additional fragmentation scale to the problem. In- 
stead, we apply the prescription suggested in Ref. [3T] 
(see also Ref. [35]), demanding 



1 — cos R 



i,Ri-,<R 



1 — COS Sq 



Pr 



VIKSn 



(2) 



where the index i runs over all partons in a cone around 
the photon of size R. For the cut-off parameter, which 
determines the QCD-IR-safe cone size around the pho- 
ton, we choose 6q — 0.7. This is a rather large isola- 
tion compared to the experimental resolution capabilities 
(e.g. Ref. [33]). Hence, the phenomenological impact of 
the full jet-photon fragmentation is expected to be small, 
in accordance with the results of Refs. [7] IMl 155] . 

We cluster all final state partons with \y p \ < 5 to jets 
via the kx algorithm [36] using a resolution parameter 
D = 0.8, adding the four momenta of clustered partons. 
The jets are required to lie in the rapidity range \yj \ < 4.5 
with transverse momenta pip > 20 GeV. The photon and 
the charged lepton are required to be hard and central, 
p e T > 20 GeV (10 GeV at the Tevatron), > 20 GeV 
(10 GeV at the Tevatron), |f#|,|?7 7 | < 2.5, while being 
separated in the azimuthal angle-pseudorapidity plane by 
Ri y = (A<^ 7 + At7 2 7 ) 1/2 > 0.4. For the separation of the 
charged lepton from observable jets, we choose Ry > 0.4 
and we require i? 77 > 0.4 for the diphoton separation. 



4 




FIG. 4: Differential maxp 3 T distribution for W^~77+jet pro- 
duction. The chosen cuts and scales are described in text. 
The horizontal line in the lower panel displays the K factor 
for total inclusive production, Tab. |T] 



FIG. 6: Differential diphoton R separation distribution for 
W" r_ 77+jet production. The chosen cuts and scales are de- 
scribed in text. The horizontal line in the lower panel displays 
the K factor for total inclusive production, Tab. [I] 
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FIG. 5: W77 invariant mass distribution for W~ -yy+jet pro- 
duction. The chosen cuts and scales are described in text. 
The horizontal line in the lower panel displays the K factor 
for total inclusive production, Tab. [T] 



Besides the photon-parton isolation criterion mentioned 
before, we also require a separation between photons and 
identified jets of R y j > 0.7. The cross sections and total 
K factors for a dynamical scale choice fiR = fip = mw-yy 
are shown in Tab. [i] Here, mw-yy denotes the W ± ^ 
invariant mass. 

The production cross section at the Tevatron is too 
small to be of evident phenomenological importance 
when viewed in the light of a total accumulated data 
set of ~ 10 fb _1 per experiment. 

We compute total K factors of 1.43 (1-48) for 
VF+77+jet (I4 y ~77+jet) production at the LHC. These 
values are quite typical for multiboson+jet production 
as found in Refs. [5-8J. The scale dependence of the 
I4 y+ 77j and W^jjj production cross sections turn out 
to be modest: when comparing /I/? = [ip = ^mwyy 
for £ = 0.5 and £ = 2, we find differences of 10.8% 
(12.0%), respectively. The dependence is dominated by 
the renormalization scale, as can be seen in Fig. I] This 
is a consequence of additional jet radiation being impor- 
tant for the probed small gluon momentum fractions. 
The right panel of Fig. [3] also shows the effect of ve- 
toing a second hard jet. The real contributions become 
smaller and there is a cut-dependent partial and acciden- 
tal cancellation within the different contributions, which 
should not to be taken as a stabilization of perturba- 
tion theory, however (see also [3 18] and below) . For the 
numerical evaluation, we split the virtual contributions 
into fermionic loops (Virtual-fermionbox, corresponding 
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to the diagrams sketched in Fig. [2]) and bosonic contribu- 
tions with one, two and three electroweak vector bosons 
attached to the quark line, i.e. Virtual-box, Figs. [T] 
(a) and (b), Virtual-pentagons, Fig. [I] (c) or Virtual- 
hexagons, Fig. [lj (d), respectively. This procedure al- 
lows us to drastically reduce the time spent in evaluating 
the part containing hexagon diagrams, which requires the 
largest amount of CPU time. The bosonic contributions 
are not individually QED gauge-invariant. However, this 
poses no problem since for our choice of gauge (effectively, 
the Coulomb gauge in the lab-frame is used for external 
photons) there are no sizable cancellations among the 
different contributions. The fact that in Fig. [3] the differ- 
ent virtual contributions share the same scale dependence 
corroborates our gauge choice. 

The phase space dependence of the QCD corrections is 
non-trivial and sizable, as can be inferred from Figs. 4]|6 
where we again choose [Ar = fiF = mwyy Additional 
parton emission redistributes the transverse momentum 
spectra. The leading jet becomes slightly harder at NLO, 
an effect which is best seen in the dynamical K factor (ra- 
tio of NLO to LO distribution) as shown in the bottom 
panel of Fig. [4] The effect on the VF77 invariant mass is 
even more pronounced. While this qualitative behavior 
is expected from kinematics, e.g. due to a photon pick- 
ing up the recoil from additional parton emission, the 
quantitative result is very important. An excess in the 
photon's transverse momentum or in the IF77 invariant 
mass at large values is easily misinterpreted as an effect 
of anomalous electroweak trilinear or quartic couplings 
|H1 132] arising from new interactions beyond the SM. 

The sizable impact of QCD corrections at large invari- 
ant masses is also visible in the diphoton separation of 
Fig. [6] At large values, when photons are highly sep- 
arated in pseudorapidity, the dynamical iGfactor rises 
well above the average value of 1.43. This experimen- 
tally clean and well-reconstructable configuration typi- 
cally amounts to a large momentum transfer in the quar- 
tic and trilinear vertices in Figs. IT] (a) and (b) and there- 
fore potentially accesses new interactions at scales much 
larger than my/. 

In Figs. [4] and [5j we also plot distributions with a veto 
on a second hard jet: no such jet with pxj > 50 GeV is 
allowed. It can be observed in Fig.[I]that the vetoed con- 
tribution does not give a sensible result for large values 
of pt of the harder jet, where large logarithms involving 



PT,cut as the other relevant scale appear: changing values 
from fj,R — iip — l/2mn/ 77 to [Ir = = 2my? 11 in- 
creases the high px tail of the distribution by a factor 2.5 
or more at NLO. In fact, in Fig.jijfor msxp T ~ 100 GeV, 
the scale-varied distributions intersect, which is yet an- 
other clear indication of the previously mentioned acci- 
dental stabilization of the vetoed cross section, which is 
cut dependent. In contrast, in Fig. [5j events with high 
invariant mass can be generated, where the high mass has 
its origin purely in the leptonic sector. These events can 
have two fairly soft jets which are not cut away, which 
yields smaller variations in the iGfactor, with changes up 
to a factor 1.3 when increasing fiR = j.ip by a factor 4. 



IV. SUMMARY AND CONCLUSIONS 

We have calculated the full NLO QCD corrections 

(-) (-) 
to the processes PP-> f ^77 + jet + X and P P — >■ 

l + 1^77 + jet + X. All off-shell and finite width effects 
have been properly taken into account. This is the first 
NLO computation which falls into the three gauge boson- 
plus-jet category. 

Quite typical for the multiboson+jet production modes 
we find large total K factors of order 1.4 for inclusive 
production, which are driven by additional jet radiation 
being significant for our inclusively chosen cuts. This 
enhancement is considerably larger than naive expecta- 
tions from a LO scale variation. The corrections exhibit 
a non-trivial phase space dependence, which could easily 
be misinterpreted as non-Standard Model physics unless 
the differential QCD corrections are properly included in 
experimental analyses. 
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